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We study the localization in the Hilbert space of a modified Tomonaga-Luttinger model. For 
the standard version of this model, the states are found to be extended in the basis of Slater 
• determinants, representing the eigenstates of the non-interacting system. The linear dispersion 

which leads to the fact that these eigenstates are extended in the modified model is replaced by one 
with random level spacings modeling the complicated one-particle spectra of realistic models. The 
£SJ ■ localization properties of the eigenstates are studied. The interactions are simplified and an effective 

one-dimensional Lloyd model is obtained. The effects of many-body energy correlations are studied 
numerically. The eigenstates of the system are found to be localized in Fock space for any strength 
of the interactions, but the localization is not exponential. 
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I. INTRODUCTION 



The concept of Anderson localization in the Fock space of many-body systems has been recently adopted for the 
studies of interacting electrons in finite-size conductors. The traditional field of localization consists of the investigation 
of electronic motion in disordered solids. These studies were initiated by Anderson jjj who considered tight-binding 
Hamiltonians on real-space lattices with random on-site energies. The behavior of the system was found to depend 
on the parameter Z = W/V, where W is the typical variance of the random potential and V is the hopping matrix 
element. For Z above some critical value the eigenfunctions of the Hamiltonian are exponentially localized with a 
characteristic scale £, called the localization length, which leads to the absence of diffusion of electrons. The Anderson 
problem inspired numerous studies (for reviews see The scaling theory of localization || predicted the strong 

dependence of the localization properties on the dimension of the lattice. It was found, for instance, that in one 
spatial dimension the states are always localized for any strength of the hopping. 

The ideas of localization can be extended to the investigation of the properties of interacting many-body Hamilto- 
nians in the basis of eigenstates of the corresponding free system. The central problem considered there is the mixing 
of a particular state representing one particle excited above the Fermi surface with many-particle states responsible 
for emergence of a quasi-particle described by the usual Landau Fermi liquid theory. For this purpose the analogy 
between the one-particle and the many-body problems was demonstrated in |jj by mapping the Hilbert space of inter- 
acting electrons onto an effective tight-binding model. Each site of the effective lattice represents the eigenstate of the 
non-interacting system and the bonds represent the interaction matrix elements. After several approximations, which 
consist mainly in neglecting the possibility of closed loops on this effective lattice, the resulting model was identified 
with the Cayley tree with branching number K depending on the energy e of the initial excitation, representing a 
particle excited above the Fermi sea. The Cayley tree model was previously studied and was found to undergo 
a localization transition for the value of parameter Z = Z c « K hxK. In the context of the original many-electron 
problem, the system is in the localized regime of the many-particle states on the states of the non-interacting system 
for energy below the lower threshold e** = Ay/g/hig, where A is the one- particle level-spacing and g is dimensionless 
O ■ conductance, i.e. conductance in units of e 2 /%. It was also argued that until the excitation energy e reaches another 
critical value e* = A^fg > e** the states are not completely delocalized, i.e. they do not mix large number of non- 
• ?-h , interacting states, but rather a small portion of them, corresponding to a subtree of the whole lattice. The last result 
was not confirmed in the subsequent studies ||] of this model by super-symmetry methods. It was found there that 
the states are strongly correlated above the threshold e > e**. 

In more recent numerical studies of the interacting models |],[l(J it was observed that the choice of the type of 
the lattice affects the localization-delocalization transition point. The crossover between the spectral statistics and 
the properties of eigenstates such as inverse participation ratio were shown to be sensitive to the type of the lattice. 
It was doubted that the actual transition and the properties of states in the delocalized phase can be accounted for 
properly using the pure tree-like structure of the lattice. It was claimed in |ll[] that the infinite tree-like lattice used 
in |fj] may be inadequate for description of the splitting of the quasi-particle peak into many-body constituents. The 
fmiteness of the lattice used in this work results in the absence of the sharp transition from localized to delocalized 
states. In the numerical work jlj], where the possibility of closed loops on the lattice was taken into account using the 
random-matrix theory no delocalization transition was observed. Rather the participation ratio and spectral width 
of the eigenstates of the interacting system exhibited a smooth crossover from almost localized to delocalized states. 
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The exact solution for the delocalized phase (which corresponds to the strongly interacting regime) is not available 
in these studies in any limit . It would be useful to have a model which can clarify the nature of the extended 
states in the strongly interacting regime. The candidates for such a system are of course one-dimensional interacting 
system whose analytical solution is available for any strength of the interactions. The present work was inspired 
by the observation that the eigenstates of the simplest interacting problem, namely the chiral Tomonaga-Luttinger 
model 13 are extended in the effective lattice of the non- interacting Slater determinants. This model is usually 
solved by the bosonization technique, which can be viewed as a generalization of the solutions by Fourier transform 
of tight-binding models without disorder. This analogy can be traced more explicitly if one considers the Hilbert 
space spanned by the Slater determinants as an effective lattice where each site is labeled by partition of some 
integer. Due to the conservation of total momentum the Hamiltonian of the chiral Tomonaga-Luttinger model is 
block-diagonal in the total momentum index N and the linear dispersion implies that energies of all the states inside 
one block are degenerate and equal NA. The interactions lift this degeneracy mixing the Slater determinants into 
coherent superpositions which are described by the bosonic quantum numbers. Thus the bosonic eigenstates of the 
full interacting model can be viewed as a generalization of the extended Bloch states in the tight-binding models with 
degenerate on-site energies. 

The linearization of the dispersion is the approximation which allows to solve the problem of interacting electrons 
in one dimension. The deviations from the linearity mixes the bosonic states, so they cannot any longer provide 
the solution to the model. In this work we introduce the random dispersion into the Tomonaga-Luttinger model 
and study the properties of its eigenstates. It is expected to model a system where the single particle spectrum is 
complicated. The question we address is whether this model exhibits localization, since the randomness introduced 
into the dispersion leads to randomness of the diagonal matrix elements of the Hamiltonian analogous to the diagonal 
disorder introduced into tight-binding Anderson models. We have found that the states remain localized in the basis 
of Slater determinants, i.e. resemble the states of the non-interacting systems. 

The outline of this paper is as follows. In section 2 we review some known facts about the standard Tomonaga- 
Luttinger model. We describe the structure of the Hilbert space as an effective lattice and discuss the bosonization 
solution in this context. In section 3 we introduce randomness into the dispersion law of the Tomonaga-Luttinger 
model and simplify the interaction matrix. Then we solve the resulting model with assumption of independence of 
the on-site energies. Section 4 is devoted to the study of relevance of the correlations between the energies. In the 
last section the results are discussed and further investigations are proposed. 



II. HILBERT SPACE OF TOMONAGA-LUTTINGER MODEL 



In this section we review a model exactly solvable for any interaction strength. The simplest model of interacting 
fermions is the one-branch (chiral) Tomonaga-Luttinger model p3| . Let n label the one-particle orbital with convention 
that in the ground state |0) of the non-interacting model, n — 1 denotes the first unoccupied state. Since a very 
large number of states is occupied (the number of fermions is large) the number n can be allowed to run from — oo 
to +oo. Another consequence is that the one-particle spectrum is assumed to be equidistant: E n = nA. Both the 
linear dispersion and the fact that the spectrum is unbound are crucial for the solubility of the model. Introducing 
creation and annihilation operators cjj, c„ for a fermion on the orbital n, with usual anti-commutation relations 
{c n , c n /} = 3n,n' (other anti-commutators vanish) the Hamiltonian of the problem can be written as 

H = H Q + U = E + J2 E n ■ 4cn - + V ™ C i+m C ncl„ m Cn> (1) 

n— — oo m— 1 n,n' — — oo 

The normal ordering of operators having infinite expectation values is defined as : O := O — (0|O|0). As a result the 
constant Eq absorbs the (infinite) free ground state energy together with the Hartree term: 

E Q = (0\H\0) = ]T £„<0| C t c „|0) + y E (0|4<v£'<vl°> (2) 

n<0 n,n'<0 

The operator of "total momentum" 

+00 

P= ^2 n : c) n c n : (3) 
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commutes with the full Hamiltonian. Due to the linear dispersion this operator is proportional to the free part of the 
Hamiltonian, i.e. the first term on the RHS of ([[]). 

Instead of solving this model by the standard method of bosonization using identities among second quantized opera- 
tors (see for example |]l3| ) we prefer to reformulate our interacting problem in order to relate it to the problems usually 
discussed within the framework of Anderson localization. We start by labeling the eigenstates of non-interacting prob- 
lem, which are the familiar Slater determinants, by shifts of particles with respect to the non-interacting ground state 
as shown in the Fig. [I} the uppermost particle is shifted up Ai levels, the second one is shifted up A2 levels etc. The 
numbers Ai are defined to satisfy Ai > Aj+i for each i. If a particle is not shifted it is not recorded. The set (A) 
of positive integers satisfying Ai = iV forms the set of partitions of N and is conveniently displayed by a Young 
diagram — array with Ai boxes in the i-th row. Following the standard notations in the mathematical literature the 
numbers Ai will be referred to as parts in what follows. For example, the partition of N = 28, (A) = (7,6,5,5,3,2) 
sometimes written as (7, 6, 5 2 , 3, 2) that is displayed in Fig. || is represented graphically as 



w - 

This pictorial representation of Slater determinants has nice properties under particle-hole interchange: recording 
shifts of the holes (which are moved downwards in energy) rather than particles in Fig. [I] yields the partition (A) = 
(6, 6, 5, 4, 4, 2, 1) whose graphical representation is just the transpose of the Young diagram for (A): 



(A)- 

Partition (A) is said to be conjugate or dual to (A). This duality results from the fact that each move of a particle 
upwards is accompanied by a move of a hole downwards. 

The eigenstates of the non-interacting Hamiltonian labeled by partitions of TV are degenerate and their energy is 
given by 



(Ei-i+xt - Ek-t) = A ^ A, = NA 



(4) 



The interactions conserve this number. The Hamiltonian (|l|) is therefore block-diagonal and does not mix the partitions 
of different N so we can restrict ourselves to the subspace of the Hilbert space spanned by partitions of N. 

The dimensionality of this subspace is given by p(N) — the number of partitions of N. There is no closed analytical 
expression for this number but the Euler expression mM for the generating function 



OO 

E 

JV=1 



p(N)q 



N 



00 

1=1 



(5) 



allows to obtain the asymptotic behavior of p(N) for large N and it is given by the formula of Hardy- Ramunajan: 



p(N) 



4V37V 



(6) 



It demonstrates that the number of "sites" in the subspace increases rapidly with N. 

Since for a given subspace N the free part of the Hamiltonian is proportional to the unit matrix the solution consists 
in diagonalization of the interacting part U. The interaction term in ([!]) is conveniently represented in our picture 
as effective binding between the sites. In order to present the results in a compact and transparent way we drop the 
requirement on a partition labeling the basis state to be represented by a non-increasing sequence. Instead we extend 
the notion to any finite sequence of integers (either positive or negative or zero) and use the following rules to relate 
this sequence to the standard (non-increasing) partition: 
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1. In any sequence (Ai, A2, . . •) two consecutive parts may be interchanged provided that the preceding part is de- 
creased by unity and the succeeding part increased by unity, the state vector which corresponds to this partition 
acquires a minus sign, i.e. 

|(Ai, .. .,Xi, Aj+i, ...)) = — I (Ai, . . ., Ai+i - 1, A, + 1, . . .)) (7) 

2. // any part exceeds by unity the preceding part the partition corresponds to zero vector, i.e. 

|(A 1 ,...,A i ,A i + l,...)) =0 (8) 

3. If the last part of the sequence is negative the partition corresponds to the zero vector. 

These rules, used in the theory of symmetric functions jl5f are in fact nothing but a standard properties of Slater 
determinants under action of second quantized operators. Rule 1 results from the fact that i-th particle is excited to 
the level rii such that 

rii = X{ + 1 — i 

with the energy Ei = n^A while the (i + l)-th particle is excited to the level 

Hj+i = Aj+i — i 

and interchanging consecutive rows in a (Slater) determinant results in a sign change. Rule 2 results from the fact 
that Ai+i = X{ + 1 implies n^x — rii which contradicts the Pauli principle. Rule 3 reflects the fact that negative Ai 
describe (forbidden) transitions to occupied states inside the Fermi sea. 

Let us denote by T)/\-\ the factors — 1, or +1 acquired when bringing a partition (A) to its standard form. In these 
notations the matrix element of the interactions between the partitions (A) and (/z) can be represented as 

U {X)M =<(A)|#|(M)> =E ^mp-Aj-m. ..)|(. .. N -m...)). (9) 

m i,j 

Since the basis of Slater determinants is orthonormal the following expression is obtained 

U WM = X/Kn2^»7(...A i -TO...) 7 ?(---M3-m...)^(...A j -m...),(...Mj-m...)- (10) 
m i,j 

The eigenstates of U are 

IG)> = X>(ol(*)> ( n ) 

(A) 

where Xm are the characters of the irreducible representation (A) of an element of the symmetric group Sn belonging 
to the class (I). For this purpose it is shown in Appendix A that 



V(...\i-m...) r l(...\i-m...\ 3 +m...)X(i) — mn mX(|) U^J 



(...Ai— m...\j+m...) (A) 

l l(...X i -m...)'l(...X i -m...\ j +m...)X "'" 

where n m is the number of times the cycle m appears in the class (I). This leads to 



m (A) \ TO / 



|(0). (13) 



(X) 

The state is therefore an eigenstate of the Hamiltonian (Q) with the eigenvalue 

E (l) =E + NA + J2 V m mn m . (14) 

n 1 

The orthogonality of the characters implies that the states form an orthogonal basis, namely: 
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{(l)\(l')) = J-S mn (15) 
9(1) 

where g = AH is the total number of elements in Sn and the number of elements in the class (/) is 

9(D= m f{ n 1 r ( 16 ) 

u X - L m™ m n m ! 

m— 1 

We normalize the states |(Z)) as 

ko> = > . \:;Vmai: (i7) 




so that they form an orthonormal basis. These are standard bosonic states usually labeled by the occupation numbers 
n m that are the eigenstates of the interacting Hamiltonian (Q) . The transformation with the matrix 

(wm = J^fx^ as) 

is the unitary transformation of the basis vectors in our Hilbert space, which diagonalizes the interacting Hamiltonian 
(|l|). When the free part of the Hamiltonian is degenerate in the subspace N the states play the role of the extended 
states (in the Hilbert space). In order to verify this statement we calculated numerically the inverse participation 
ratio P, defined as 

(A) 

for some arbitrary states (/). Although no closed expression is available for the general amplitude (|lj) of the state 
on the site (A), the beautiful recurrent relation of the characters of the symmetric group ( A12| ) (based on the Frobenius 



formula) was implemented numerically to calculate these amplitudes. The inverse participation ratio averaged over 
the bosonic states (I) is shown in Fig. |||for different values of N. As the number of sites N S i te — p(N) increases the 
inverse participation ratio P was found to decay as P = 1/Nf ite , where a w 0.8. This allows us to conclude that the 
bosonic states (I) are indeed extended in the basis of free Slater determinants (A). 



III. SIMPLIFIED INTERACTIONS 



The existence of extended states in the chiral Tomonaga-Luttinger model can be attributed to the degeneracy of 
the diagonal elements of the Hamiltonian (|l|), since the interactions ( JlOj ) have short-range in the space of partitions. 
In order to define more precisely what we mean by short-range in our Hilbert space we have to order the partitions 
in some order. Then the interaction can be shown to connect only blocks of elements that are not too far one from 
the other. 

It is natural to order the partitions in the lexicographic order, i.e. if = 1,2,3... numbers the letters in an 
alphabet the lexicographic order of partitions corresponds to the order of words in Arabic or Hebrew dictionary. For 
example, the 11 partitions of A^ = 6 in the lexicographic order are: 



k = 


1 


(1,1,1,1,1,1) 


k = 


2 


(2,1,1,1,1) 


k = 


3 


(3,1,1,1) 


k = 


4 


(2,2,1,1) 






(4,1,1) 


k = 


5 


(3,2,1) 






(5,1) 


k = 


G 


(2,2,2) 






(4,2) 






(3,3) 






(6) 



(20) 
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The partitions are grouped by fc — the number of l's (that appear always on the right). The fc-th group contains 
partitions of the form 

(A 1 ,A 2 ,...A i ,l Ar - fc ) (21) 

with constraint Xi > 2, with the only exception for the first group containing (1^). The typical Young diagram of a 
partition in the fc-th group is of the form 



N-ki 



(22) 



To calculate the number of partitions in the k-th group we note that exactly p(k) partitions have at least N — k ones. 
Among them p(k — 1) have at least N — (k — 1) ones and the rest n(k) = p(k) — p(k — 1) have exactly N — k ones. 
Therefore the fc-th group contains n(fc) partitions. 

Usually the matrix element of the interaction potential V m is a decreasing function of the momentum transfer to, 
with some effective range m c . In what follows it will be convenient to analyze a model where the only non-zero 
matrix element corresponds to m = 1 with V\ = V. It will not affect qualitatively our results, since we are primarily 
interested in the limit of large N subspace. The matrix element of this interaction is given by the simplified version 
of© 

U W,W = ^XX..A i -i...),(...^-i...)- ( 23 ) 

The factors r\ are omitted since the partitions (. . . A, — 1 . . .) and (. . . fij — 1 . . .) are either already in the standard 
(non-increasing) form or corresponds to zero vector and do not contribute to the sum. Therefore the interaction 
couples partitions (A) and (ji) that differ in two of the Aj's, namely Xi — fii = 1 and Xj — fij = — 1. The Young diagram 
of the partition (A) is obtained from the one of the partition (ji) by moving one square under the constraint that the 
new Young diagram is legitimate. Consider the fc-th group. If fii > 1 and fij > 2 the coupling is between partitions 
inside this group. This corresponds to moving a square between \ii > 1 rows of the Young diagram. In order to couple 
to neighboring groups a square corresponding to fii — 1 should be moved. Moving a square from/to any row with 
/ij > 1 to/from the region where fii = 1 leads to a coupling to the neighboring group k ± 1. The exception is if a 
square is moved from fii = 2 (that has only two squares and after the move only one is left) and is attached as the 
last square in the tail of fij — I, or if the last square is moved to fij = 1 turning it into a fi 3 ■ = 2 row. In this case the 
length of the tail of fij — 1 is changed by two, i.e. the coupling is to the group fc ± 2. These are the only moves of one 
square, therefore only the groups fc' such that \k — k'\ < 2 are coupled to the fc-th group and the interaction between 
the various blocks of partitions is of short range. Fig. || shows the structure of the interaction matrix for N = 14. 

It is therefore natural to expect that the diagonal disorder will destroy the extended states of the chiral Tomonaga- 
Luttinger model. In the original model in the absence of interactions the eigenenergies E^x) are independent of the 
partition and take the constant value E^x) — Eq + NA. We introduce a modification and assume that the E(X) are 
independent random variables. This assumption is made for the sake of simplicity. It cannot be correct since there 
are p(N) values of E(x) that depend on 2N level spacings. An assumption that is more physical will be introduced in 
the following section. 

In order to be able to solve the model with random energies we simplify further the interaction (|23|). From the 
experience in research of Anderson localization we expect that the behavior of the eigenvectors will not be affected 
by the precise form of the short-range off-diagonal matrix elements, so we make the following simplifications. We 
assume that the effect of transitions inside the same group of partitions is to delocalize the states within the group. 
The localization within a group is expected only to enhance the overall localization. Next, we neglect the coupling to 
the fc + 2-th group. The reason is that for a given partition (A) there exists only one coupling term of this kind. It 
couples 

((...), l N ' k ) and ((...)', 2, l N - k - 2 ) (24) 



G 



This should be compared it with the couplings of (A) to the next group: 



((...),i JV - fc )^((-.-)M JV - fc - 1 ) 



(25) 



The number of such couplings is given by the number of different A^ > 2 in the partition (A), which is generally much 
larger than 1 for large enough k. It is equal to the number of rows with Ai > 2 in the Young diagram of the partition. 
Consider for example the partition (6, 4 2 , 3, l 4 ) of N = 21 belonging to he the 17-th group corresponding to the Young 
diagram (§2h. It has (apart of ones) three different parts: 6,4,3 and is coupled to 3 partitions in the 18-th group: 



(26) 



(7,4 2 ,3,1 5 



(6,5,4,3, l 3 ) (6,4 3 ,1 3 ) 



Another way to visualize the number of coupling is to count the number of concave corners in the Young diagram 
that is equal to the number of parts namely lines of different length larger than unity. 

In Appendix B we show that d{k) — the number of different parts in all partitions of k with A^ > 2 is exactly 
p(k — 2) — the number of (any) partitions of (k — 2). Therefore the mean number of transitions from the fc-th group 
to the k + 1-th group is given by t(k) — d(k)/n(k) — p(k — 2)/ (p(k) — p(k — 1)). 

The Schrodinger equation for our simplified model can be written as a tight-binding form where the groups of 
partitions are the effective sites: 



O) 



(27) 



where ip{\) k denotes the value of the wave function on site(A) in the fc-th group. Let us assume that the wave function 
is extended within the kth group. We introduce a simpler model that is expected to exhibit similar localization 
properties fl7| . We replace the interaction ?7(a).(/j) by an average interaction that couples each partition in the fc-th 
group to all partitions in the group k ± 1, 



n(k + 1) 



In this approximation (|27|) reduces to 



(E - E w )1>(\) k = V k ^ V ip(fJ,) k -x +V k Y] i>{v) 



(/*) 



>k+l 



(28) 



(29) 



The numerator in (|2q ) can be estimated as V times the mean number of couplings with partitions in (k + l)-th group 

t(k) _ T/ d(k) 



V k = V- 



= V- 



(30) 



n(k + l) n(k)n(k + l) 

Denoting the sum of the wave- functions over a group by C k — E(A) we nn< ^ w ith the help of (^9|) that 

£kC k — Vk-iCk-i + VfcCfc+i (31) 

where 

£ k 4^ E-E, 



(A) 



(A) 



(32) 



and the sum is over all partitions in the group k. 

For random independent all the £ k are independent random variables and the tight-binding model ( pl| ) is an 
Anderson model at zero energy. In the spirit of transfer matrix method it is expected that also for £ k given by (p52) 
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similar localization properties will be found. The properties of the solution of this one-dimensional model can be 
analyzed within the framework of Anderson localization in one-dimensional systems [fjfjf . 

One can take the uniform distribution of energies that turns out to be convenient for calculations. Remembering 
that each many-body energy is a sum of N level spacings we assume the following distribution: 



P(E) = { VV^A, -VNA/2 < E < VNA/2 
1 0, otherwise 



(33) 



For a variety of distributions the distribution of £k tends to the Cauchy distribution: 

1 h(E) 



p < £ <> " (34) 



where 8 k (E) = 8{E)/n{k) and for (g3[) 6(E) = VNA/tt. 

Multiplying the equation (|3|) by n(k)n(k + \)/{Vd{k)) and observing that in the limit k — > oo one finds n(k)/n(k + 
1) — ► 1 and d(k)/d(k — 1) — > 1 so that (|3l|) takes the form 

ffcCfc = Cfe-i + Cfc+i (35) 
It is a Lloyd model with a fc-dependent distribution of on-site energies characterized by parameter 

6(E) n(k+l) 



which in the limit k — > oo has the property 



[N A 



We generalize the solution of the Lloyd model [M to the present case of site-dependent distribution of energies. We 
observe that the distribution parameter depends weakly on the site number k and use the site-dependent localization 
length ansatz: 



i/& = m ^ = cosh- y i + f « m 4 = m y - - (38) 

The last approximation is justified in the weak-coupling limit V/A <C 1. 

The validity of the site-dependent localization length assumption is checked numerically. The on-site energies £k 
were generated from the Cauchy distribution (|34|) with parameter 6k given by (|37]). The inverse localization length 
l/£ (Lyapunov exponent) was calculated using the standard technique of transfer matrices [pof and it was found that 
the site-dependent ansatz (^) describes well (at least qualitatively) the behavior of the localization length. Thus for 
independently distributed energies the typical behavior of the wave functions is described by coefficients that decay 
as 

C k ~ e- k '^ = exp (-k log ^/J^ ■ (39) 
In the strong coupling limit similar considerations lead to 



but in this limit these require some further justification. The importance of this result is that numerical studies of 
more realistic many- particle energy distributions confirm the validity of ([59]). This will be the subject of the next 
section. 
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IV. EFFECTS OF ENERGY CORRELATIONS 



In the last section we considered a model in which the diagonal matrix elements -27(a) are independent random 
variables. We did not specify the distribution law of these energies, since the solution does not depend on fine 
details, but rather on general properties of the distribution. We would like to stress that the assumption of statistical 
independence of the energies Bai is not realistic for a many-body system. To illustrate this statement a simple 
argument can be given. Suppose we are dealing with a disordered mesoscopic system and the one-particle spectrum 
is random, i.e. consists of levels with independently distributed spacings Aj — Ej — Ej-\. The generalization of the 
expression (||) for the energy of the partition (A) for this case is 

Xi+l-i 

E w = Y, ( E i-i+\-Ex-i)=Yl E A ^ ( 41 ) 

i i j=2— i 

In order to describe the energies E^x) in the subspace N we need 2N random quantities Aj. On the other hand 
the number of states in this subspace is given by p(N), the number of partitions of N which grows exponentially for 
N — > oo. Therefore there exist much more many-body energies E^x) than independent level spacings. Therefore the 
many-particle energies E(x)i which are the sums of the one-particle energies, are statistically dependent. 

The analytical calculation of the (joint) distribution of the parameters £ k defined in ( |32| ) for the effective tight- 
binding model ( pi] ) seems to be hopeless. Nevertheless we expect some features of this distribution, based on the 
formula (|32|). For an eigenstate the eigenvalue E is in the interval [—1 + EV^\ n , +1 + E7j? x ], where EJ^ n and E7j? x 
are the minimal and maximal values of -E(a) ■ The value of E can be very close one of the Erx) ■ Therefore the terms 
in the sum on the RHS of ( |32|) can be positive or negative and they strongly fluctuate in magnitude. Therefore there 
is a finite probability that this sum will take a value in the vicinity of zero. Consequently it is expected that the 
distribution P(£ k ) is broad with a diverging second moment or variance: 

A£l = ((£ k - (£ k )) 2 ) -> oo (42) 

Another important issue is the effective statistical independence of £ k for different k-s. It is known for Anderson 
localization that the long-range correlations between the on-site energies change dramatically the criteria for the onset 
of the localization (see pq | and references therein). The measure of statistical dependence of energies in our case is 
the correlation matrix defined as 

((£ k - (&))(£, - (£,))) 
° kl = A^Ai, (43) 



Because of the strong fluctuations between the terms of the sum on the RHS of (|32j), that arc different for the various 
groups, it is reasonable that the off-diagonal terms of the correlation function are much smaller than the diagonal 
ones. Usually for localization only pair correlations are important [ fl9|| . In the view of the diverging variance we 
assume that C k i = Ski > an assumption that will be tested numerically. 

These facts were checked numerically and in the approximation of independent energies the effective distribution 
P(£ k ) can be calculated. For this purpose M = 1000 realizations of N = 30 level spacings Aj were generated from 
the exponential distribution: 



P(A) = e- A (44) 
30) = 



with A = 1. Then we generated numerically the whole list of jp(30) = 5604 partitions of N — 30 and for each 
realization of level spacings the energies E^x) were calculated using (hTj). Out of these energies the effective energies 
£ k were obtained using the definition (j3^) for E = 0. Several statistical tests were performed and results are shown 
inFig.||. 

The correlation matrix of the energies £ k was computed. The (. . .) averages were performed over realizations of the 
level spacings. The result is shown in Fig. ||. It is plausible that the correlations are of short range in the space of 
indices k and only the diagonal elements of C k i are of appreciable magnitude. This behavior can be attributed to the 
diverging variance. 

Due to the fact that the distribution of the energies £ k is expected to have a diverging variance the distribution 
P(£ k ) cannot be calculated from a histogram. Instead its characteristic function 

X k (p) = (e ip£ *> (45) 







was computed and plotted for several values of k in Fig. || as a function of p. The diverging second moment due to 
the fact that P(x) oc 1/x 2 manifests itself in the discontinuity of the derivative at p = as can be seen from the Fig||. 
Therefore we have strong evidence that the energies £ k have a broad distribution. 

Observing the behavior of the characteristic function shown on the Fig. |H|and Fig. |^ we introduced a scaling ansatz 
for the distribution function: 

P(S k )=P k (£) = ^p(n(k)^y (46) 

The universal function P(x) is independent of k. The scaling of the distribution implies that the characteristic function 
( pJ5| ) satisfies 

X k (p) = X , X(q) = r dxP(xY qX (47) 



n(k) 

which is easy to check by plotting X k {p) as a function of q — pA/n(k) for different values of k. These curves are 
displayed in Fig. |7] and are indeed found to be close to each other in spite of the fact that n(k) ranges from n(l) = 1 
to n(30) = 1039. 

We found numerically that the universal characteristic function can be well approximated by the general form 

X(q) = cxp f-a\q\ - b 4-) (48) 



The value of parameters were found to be a ~ 1.15, b ~ 6.32. This is the characteristic function of a distribution of a 
sum of two random variables x\ and x 2l where x\ is drawn from Cauchy distribution with probability density 

7T x A + a z 

and the second is distributed normally according to 

1 _ai 

VZTTO 

The distribution of x = x\ + x-i is thus given by 

/+oo 
dypi(x - y)p 2 (y) (51) 
- CO 

and its characteristic function is pjjj). It has a diverging variance due to the fact that P(x) behaves as for large 

The scaling factor A/n(k) and the limiting behavior of the distribution function are the same as those found in 
the case of mutually independent energies. This allows us to conclude that the behavior of the coefficients C k given 
by (|39|), that was obtained in the model with uncorrelated -E'(a)) is satisfied also by the eigenf unctions of the more 
realistic model of the present section. 



V. DISCUSSION 



The localization properties of the eigenstates of the modified interacting chiral Tomonaga-Luttinger model were 
studied. For total energy NA subspace we defined an effective lattice of Slater determinants labeled by partitions of 
N which represent the eigenstates of the system in the absence of interactions. We have observed that in the case 
of linear dispersion the eigenstates of the interacting model described by the bosons are extended over this lattice. 
When the random one-particle levels are introduced instead of the equidistant spectrum, the bosons are no more the 
eigenstates and the problem of diagonalizing of the Hamiltonian. In order to achieve this goal we have grouped the 
different sites (partitions) according to the number of ones, representing the distance of the most excited particle 
from the Fermi surface in the original Slater determinant. This is to be contrasted with the grouping of many-body 
states into generations according to the number of excited particles and holes as was firstly done in |(| . Each resulting 
group of partitions was given a label k. The interactions were shown to be of short range in k. The interactions 
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were further simplified by assuming that they interconnect only groups of partitions corresponding to the adjacent 
values of the index k (and the interactions between groups k and k ± 2 were neglected). Assuming further that the 
states are extended within each group k we obtained an effective tight-binding Anderson model where each group k 
has an effective random energy and the hopping takes place between the nearest neighbors only. In the assumption 
of independently distributed on-site energies the model is found to be close to the Lloyd model characterized by a 
broad distribution of the on-site energies. The main feature of our effective model is that the distribution of energies 
is fc-dependent. We employed the site-dependent localization length assumption, which consists in assuming the 
parameters of distribution change slowly, so locally the localization length is given by the same formula as for the 
Lloyd model with parameters of the local distribution instead of the uniform global one. Up to an unknown numerical 
factor this solution describes well the behavior of the localization length as a function of k and it enabled us to show 
that the eigenstates are localized in the space of groups of partitions, k. The weak-coupling and the strong-coupling 
regimes can be identified according to the value of the parameter V/A. In the weak-coupling regime V *C A, the 

amplitudes decay as exp ^— k log \J^y^j , namely nearly exponentially. In the strong coupling regime V 3> A, this 



behavior changes to exp \ ~\J ^WV^j ■ 



One important point which was discussed is the effect of correlations between the many-body energies of the non- 
interacting system. In the studies devoted to the localization in many-body systems the effective on-site energies 
corresponding to the eigenstates of the non-interacting system are assumed to be independent random variables, so 
their spectrum is described by a Poissonian sequence. This is not the case for the true many-body fermionic system, 
since these energies are sums of independent energies of one-particle excitations. The number of independent random 
variables is therefore much less than the number of different many-body states. This situation is similar in some sense 
to one faced in random walk problems, where the increments are distributed independently, while there is a strong 
correlation between the positions of the random walk at different times. In our work we considered this problem 
numerically and found that although the correlation between the exact energies corresponding to the eigenstates of 
the free model are not a priori negligible the random energies of the effective sites k are found to be very weakly 
correlated. Moreover it was found that the parameters of the distribution of the effective energy of the group k scale 
with the number of the partitions inside the group in the same way as they do in the case of independent many-body 
energies. This led us to the conclusion that our result for the localization of the eigenstates are valid for more realistic 
assumptions on the distribution of the many-body energies. 

The main open problem arising from the present discussion is the study of the interacting electrons in dimensions 
higher than one by the methods of high-dimensional bosonization. In this technique [p4|-pq] the fermionic system is 
described as an infinite set of one-dimensional Luttinger models attached to a point on the Fermi surface. Therefore the 
interactions mix the states belonging to different points in addition to the coupling within the same point (represented 
by a Luttinger model), which can be treated within the framework developed in the present paper. This possibly can 
be a mechanism for the derealization of the states and a resulting localization transition must be compared with the 
transition encountered in the recent studies [|6|j^-|ic|l of interacting electrons in finite systems. 

Another question to be explored in the future is the nature of localization in the Hilbert space for single-particle 
dispersions that are more realistic than the random ones. 
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APPENDIX A: A PROPERTY FOR CHARACTERS OF S N 

In this appendix we show that the property ( p"2| ) 

E(...Ai-m...A,-+m...) (A) 
%..\ i -m...)V(...\i-m...\ j +m...)X\i) = ™mX( ( ) 

i,j 

is satisfied for the characters of the symmetric group Sn. We did not find this property in any textbook on the 
representations of groups, but its validity follows from the known properties of the characters of Sn found, for 
example, in jlj| and . We present the proof of this property here for the sake of completeness of the discussion. 
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The Frobenius formula as it is written in the classic textbooks on group theory 0,|T5| states the following identity 
between two antisymmetric functions of TV-dimensional vector z = {z\, Z2, ■ ■ ■ , Zn): 



S (l) (z)D(z) = J2x$*(x)(z) 
(A) 



(Al) 



where Xm is the character of the class (I) of the symmetric group Sn in the irreducible representation (A). With the 
partition (A) one associates the totally antisymmetric function: 



^ {x) (z) = ^2sgn(P)z^t ) N - 1 



-1 A2+W-2 

Z P(2) 



~P(N) 



(A2) 



where the sum runs over all the permutations of 1,2,..., TV, the sign of permutation being ± for an even/odd 
permutation. This is, up to a normalization constant, a Sl ater determinant constructed out of one-particle wave- 
functions z rii = Z x i+ N -^ where z = e' te . The factor D(z) in (|Al[) is the Vandermonde determinant: 



D(z) 



1 


1 • 


1 


Z\ 


Z2 


ZN 


,.2 


7 2 ■ 


7 2 


-1 


z 2 


Z N 


JV-1 


/-I . 


■ z N ~ x 


1 


z 2 


N 



(A3) 



representing the ground state of TV fermions. Let S m be a power sum of variables Zj defined as follows: 



(A4) 



The function S m (z) is totally symmetric function of the Zj-s. These functions are linearly independent as can be 
checked by calculating the Jacobian: 



dzj 



NlD(z) 



The class (I) is characterized by m, ri2, • ■ • , n m , . . . cycles of length 1,2,. 
function is defined as a product of power sums 



(A5) 



For this class the totally symmetric 



S {l) (z) = S^S™...S™ = 



(A6) 



Since for different m-s the power sums S — m are linearly independent so are the monomials Sm for different partitions 
(I). 



The Frobenius formula (Al) can be inverted with help of the orthogonality relation of the characters: 



(A7) 



(0 



where g = N\ is the order of the Sn and gm is the number of elements of Sn in the class (I) given by (|l6|). 

We turn now to develop a recurrence relation between the characters. For this purpose it is instructive to consider 
the product of the antisymmetric function ^r^\(z) and a power sum S rn (z). Since the power sum is totally symmetric 
the summands can be rearranged in the order set by each permutation and one obtains 



S m (^* M (*)=5>gn(P) (l>?wj z p\t) N ~ 1 

= 7 l(-^+rn...)^(...^+ m ...){z) 



fi 2 +N-2 
P(2) 



-P(N) 



(A8) 



where W... w +m...) — — 1> or +1 is determined accordingly to the rules of Section 2. 
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Having established the relation between the Slater determinants ^(u) for various partitions we turn to develop a 
recurrence relation between the characters. Consider class (I) with at least one m-cycle. Let us denote by (l' m ) the 
class obtained by removing one m-cycle from (I). It is then clear that 



S(i) — S m Sa/ 



c;, 



5>rU 



V'n, 



(A9) 



The partition (l' m ) defines a class in the symmetry group Sjv- m . The Frobenius formula (A.1) for this class reads 



fa') 



(A10) 



while for (/) it can be rewritten with the help of (A8) in the form 



S {l) D(z) = S m S {lL) D(z) = EX(fi 5 »** 



(m') 



£xg|3£»?(...,tJ+m..0*(...K 

0*') ; 



j+m...) 



(All) 



Comparing the coefficients of the ^(^') in ( |All| ) with those in ( Al) we find the recurrence relation of the characters 



(...Ai-m...) 



(A12) 



V (A (.... 

This relation yields the characters of Sjv in terms of the characters of the symmetric group Sn- 

In order to derive ( |l2|) take some particular m and let each partition (/) contain m as its part rim times. For each 
(A) consider the following sum over partitions: 



(A13) 



(0 



If (7) has no parts equal to m, rim = and (I) does not contribute to this sum. We want to factor out S rn (z) from 
the sum and for this purpose we can write Sn)(z) — S m (z)S^i ){z), so that the sum is over (l' m ). With the help of 
( filf ) we see that the expansion coefficients satisfy 



x 1 

9(i)mn m /g = mn m J| — 



TT — 



k—l k^m 



(A14) 



where g' — (N — m)\. The sum in ( A13) takes the form 



(A15) 



Application of the recurrence relation ( A12) to the characters in the sum together with the inverse Frobenius formula 
(|A7|) yields 



F{X){z)D{z) = S m (z)^>)(...j,- m ...)^3(ij);( l j ™"' ] S(v m) {z)D{z) 



9 (O 



= E ? /(...A s -m...)5'm(2)*(...A I -m...)(^) 

i 

The relation (Ia|) is applied to show that 



(A16) 
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V(...\ i —m...) r ](...\i~m...\ j +m...)^(...\ i -m...\ j +m...){ z ) 



1,3 



Application of the inverse Frobenius formula (A7) yields 



(...Xi— m...A,-+m...) a i \ 
'>l(...\i-ni...)V(...\i-ni...\ j + m ...)g(l)X(l) *(i)( Z J 



(A17) 



(A18) 



a (i) 



Making use of the linear independence of Sm for different partitions (I) and equating the coefficients of the Sm(z) 
with the corresponding coefficients in (A1S) results in the required relation (p^). 



APPENDIX B: NUMBER OF DISTINCT PARTS OF PARTITIONS. 



In this appendix we show that d(N) - the number of distinct parts in the partitions of N with parts Xi > 2 equals 
p(N — 2) - the number of unrestricted partitions of N — 2. Together with the generating function (Q) for the number 
of partitions let us introduce the generating function 



(Bl) 



JV=1 d=l 



for the number p(N, d) of partitions of N having entries that satisfy A,; > 2 and having exactly d distinct parts. For 
example for N = 6 and d — 2 this number is p(N,d) = 1, corresponding to the partition (4,2) as can be seen from 
the list (EOh . The number of distinct parts d(N) in all partitions of N with parts > 2 is then given by 



JY 



d(N) =^dp(N,d) 



i=l 



It is generated by the following function 



N=l 



z-rP{z,q) 
az 



The explicit expression for p[z, q) can be found easily summing directly over all the appropriate partitions: 

(A),A«>2 



(B2) 



(B3) 



(B4) 



where d(X) is the number of distinct parts in (A) and |(A)| = J^i ^i- It i s convenient to use the sets of "occupation 
numbers" {n m }, where n m - the number of times m appears in the partition (A). Noting that |(A)| = Y], Xj — 
J2m mn m anc I d(X) — J2m®( n m), where 9(x) is the Heavyside step function we rewrite (B4) as a sum over all the 
configurations {n m }, with m > 2 



{n m },m>2 



e: =2 



V°° etr. 

r/—tm = 2 v 



oc oo 

n e i 

rn — 2 n m — 



m=2 



n in 



zq" 



q" 



Differentiating this expression with respect to z at z = 1 we get 

d T-r / . zq 



d{q) 



dz 



m=2 



n i+r 



z=l 



g 2 p(g) 



(B5) 



(B6) 



where p{q) is given by (ph. The function q 2 p(q) is the generating function for p(iV — 2). The equality of the generating 
functions d(q) and q 2 p(q) implies the equality of the coefficients: d(N) = p(N — 2), i.e. the number of different parts 
in all partitions of N with entries > 2 is given by the number of unrestricted partitions of N — 2. Applying this 
argumentation to all partitions ending with N — k ones yields d(k) = p(k — 2). 
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APPENDIX C: DISTRIBUTION OF HARMONIC SUM OF INDEPENDENT RANDOM VARIABLES. 



We consider the harmonic mean X defined as 



N 



X ~ N 

71—1 



(CI) 



of N independent random variables x\, xi, ■ ■ ■ x n drawn from the common distribution with probability density P{x). 
The only requirement we impose on P(x) is that 



lim P(x) = C > 



(C2) 



Consider the random variable y — X/x. Its distribution function Q(y) can be calculated by the standard methods. 
For y — ► ±00 it has the asymptotic behavior 



Q(v) 



c 



(C3) 



and its second moment diverges. The general theory of broad distributions |21-23| predicts that in the limit N — > 00 
the probability density of the sum (CI) Y — 1/X is Lorenzian 



it (vrC) 2 + Y 2 

and consequently the probability density Pn(X)) of variable X = 1/Y is again a Lorenzian: 

p N (x) = L ] 



TT S 2 + X 2 

characterized by the parameter S related to the original distribution P(x) by the formula: 



(C4) 



(C5) 



(C6) 



In particular for the uniform distribution (p3n it takes the value S = A/n 
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FIG. 1. Correspondence between Slater determinant and partition (A) = (7,6,5,5,3,2). The dual partition 
(A) = (6, 6, 5, 4, 4, 2, 1) represents shifts of the holes. 
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FIG. 2. The inverse participation ratio ([lS)]) as function of number of sites N 3ite = p(N). The solid line represents the 
P = 1/N? ite , where a w 0.8 




FIG. 3. The matrix of the interaction feq). The blocks corresponding to different groups of partitions are shown. 
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FIG. 4. The absolute value of correlation matrix elements Cm defined in ( |43| ) of effective energies for N = 30. 
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FIG. 5. Characteristic function Xk(p) of the effective energy distribution for some values of k. 
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FIG. 6. The same as on the previous figure for a smaller range of p. 
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FIG. 7. The characteristic function Xk(p) as a function of the scaled variable q = pA/n(k) for values of k used in Figs 
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